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Abstract. The cold diffuse interstellar medium must harbor 
pockets of hot gas to produce the large observed abundances 
of molecular species, the formation of which require much 
more energy than available in the bulk of its volume. These 
hot spots have so far escaped direct detection but observations 
and modeling severely constrain their phase-space structure i.e. 
they must have a small volume filling factor (a few %), surface 
filling factors larger than unity with large fluctuations about av- 
erage and comparable velocity structure in "pencil beams" and 
"large beams". 

The dissipation of the non-thermal energy of supersonic 
turbulence occurs in bursts, either in shocks or in the regions 
of large shear at the boundary of coherent vortices. These two 
processes are susceptible to generate localized hot regions in 
the cold medium. Yet, it is of interest to determine which of 
them, if any, dominates the dissipation of turbulence in the cold 
interstellar medium. 

In this paper, we analyze the spatial and kinematic proper- 
ties of two subsets of hydrodynamical compressible turbulence: 
the regions of largest negative divergence and those of largest 
vorticity and confront them with the observational constraints. 
We find that these two subsets fulfill the constraints equally 
well. A similar analysis should be conducted in the future on 
simulations of MHD turbulence. 

Key words: ISM: evolution - ISM: kinematics and dynamics - 
ISM: molecules - ISM: structure - Turbulence 



1. Introduction 

For several decades now, molecules have been detected in the 
diffuse component of the cold neutral medium and these obser- 
vations raise several intrigui ng questions. Firstly, the large ob - 
served abundances of CH+ (|Crane et al„ 1995[ |Gredel, 1997| ) 
and HCO+ and OH ( |Lucas and Liszt, lfgg j in the cold dif- 
fuse medium imply that activation barriers and endothermic- 
ities of several thousands Kelvin be overcome. The formation 
of CH + proceeds through the endothermic reaction of C + with 
H 2 (Ai;/fc=4640 K) and that of OH via the reaction of O with 
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H 2 which has an activation energy AE/k =2980 K. In the dif- 
fuse gas, HCO + forms from CH;j\ a daughter molecule of 
CH + . None of the large observed abundances can be explained 
by standard steady-state chemistry in cold diffuse gas. Pockets 
of hot gas must therefore exist in the cold diffuse medium. 

Secondly, many other molecules like CO, CS, SO, CN, 
HCN, HNC, H 2 S, C 2 H have now been detected in absorption 
in front of extragalactic continuum sources in local or red- 
shifted gas (Lucas & Liszt 1993, 1994, 1997; Liszt & Lucas 
1994, 1995, 1996; Wiklind & Combes 1997, 1998). The lines 
of sight sample the edges of molecular clouds or the diffuse 
medium, which corresponds to gas poorly shielded from the 
ambient UV radiation field. This derives from the low excita- 
tion temperatures (often close to the temperature of the cosmic 
background) measured in these transitions. Yet, the molecular 
abundances derived are close to those of dark clouds. 

Thirdly, the spatial and velocity distribution of these 
molecules is highly elusive. Observations of CO emission and 
absorption lines toward extragalactic sources by Liszt & Lucas 



(1998) show that: (i) absorption and emission line profiles have 



comparable linewidths whereas the projected size of the sam- 
pled volumes of gas differ by more than four orders of mag- 
nitude (60 vs. 10 -3 arcsec), (ii) the excitation temperature of 
the CO molecules changes only weakly across the profiles and 
(Hi) there are very few lines of sight with no absorption line 
detected. The phase-space distribution of the CO-rich gas in 
diffuse clouds must therefore have a surface filling factor close 
to unity, and its velocity field must be as adequately sampled by 
a pencil beam than by a large beam. As mentioned by the au- 
thors, these observations suggest a one-dimensional structure 
for the molecular component, which contrasts with the profuse 
small scale structure observed in emission. 

Two mechanisms, operating at very different size scales, 
can trigger a hot chemistry in the cold diffuse medium 
and produce the observed molecular abundances: magneto- 



hydrodynamical (MHD) shocks (Elitzur and Watson, 1978 
praine and Katz, 1986 ; Draine, 1986; Pineau des Forts et al 



1986; Flower and Pineau des Forts, 1998| ) and intense coher- 
ent vortices, responsible for a non-negligible fraction of the 
viscous dis sipation of supersonic turbulence (Falgar one and 
Puget, 1995; |Falgaroneetal., 1995| ; |Joulain et al„ 1998| ). This is 
so because the dissipation of the non-thermal energy of super- 
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sonic turbulence is concentrated in shocks and in the regions 
of large shear at the boundary of coherent vortices. In both 
cases, only a few % of hot gas on any line of sight across the 
cold medium is sufficient to reproduce the observed column 
densities of molecules. This small fraction of hot gas corre- 
sponds to about 6 MHD shocks of 10 km s _1 or 1000 vortices 
with rotation velocity of 3.5 km s _1 per magnitude of gas (or 
AT tot = 1.8 x 10 21 cm' 2 ) of density n rj 30 cm -3 (Verstraete 
etal., 1999). 

The main problem met with the models of individual C 
shocks or vortices is that the predicted shift between neutral 
and ionized species is larger than observed. In this paper, we 
simply address the issue of the impact of the line of sight av- 
eraging upon the resulting line profiles in a turbulent velocity 
field. We investigate the spatial and velocity distribution of the 
regions of high vorticity or high negative divergence (shocks) 
in a simulation of compressible turbulence to test whether or 
not the space-velocity characteristics of any of these subsets 
fulfill the observational requirements. We carry our analysis on 
the numerical simulations of compressible turbulence of Porter, 
Pouquet & Woodward (1994). They are 512 3 simulations of 



midly compressible turbulence (initial rms Mach number=l). 
The time step analysed here is t = 1.2r ac where r ac is the 
acoustic crossing time. At that epoch in the simulation, many 
shocks have survived producing density contrasts as large as 
40, but the bulk of the energy at small scales is contained in the 
vortical modes. 

These simulations are hydrodynamical and do not include 
magnetic field. The impact of magnetic field upon the statistical 
properties of turbulence may be less important than foreseen. 
Recent simulations of MHD compressible turbulence (Ostriker 
et al., 1999) show that the dissipation timescale of MHD turbu- 
lence is closer than previously thought to that of hydrodynam- 
ical turbulence. Descriptions of the energy cascade in models 
of MHD turb ulence also predict that vortic es must play an im- 
portant role (Goldreich and Sridhar, 1995; Lazarian and Vish- 



niac, 1999). Magnetic field does not p revent either the intermit- 
tency of the velocity field to develop (Brandenburg et al., 1996 



Galtier et al., 1998HPolitano and Pouquet, 1998bfc Politano and 
Pouquet, 1998a; [Politano and Pouquet, 1995| ). In MHD as in 
hydrodynamic turbulence, shocks interact and generate vortex 
layers, which are Kelvin-Helmholtz unstable and eventually 
form vorticity filaments. 

Another limitation of our study lies in the fact that the 
Reynolds number of the simulations is small compared to 
that of interstellar turbulence. Recent studies of high Reynolds 
number turbulence in Helium bring the unexpected result that 
the statistical properties of the velocity field have little de- 
pend ence with the Reynol ds number, at large Reynolds num- 
bers ( labeling et al., 1996 ). We therefore believe that the anal- 
ysis presented here on hydrodynamical simulations has some 
relevance to the understanding of interstellar turbulence. 

We discuss the velocity structure of the two subsets (regions 
of high vorticity or of high negative divergence) and compare 
in each case the velocity samplings provided by pencil beams 
versus large beams (Sect. |J). In Sect. 0, we discuss the spa- 
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Fig. 1. PDFs of the values of (a) the positive and negative di- 
vergence and (b) the vorticity in the whole cube. The integrated 
numbers of points below (case a) above (case b) each threshold 
indicated by a vertical bar are 3%. 

tial distributions of the two subsets. In Sect. |], we compare the 
results derived from the numerical simulations with the obser- 
vations. 

2. Velocity structure of subsets of intense vortices and 
shocks in compressible turbulence. 

2.1. Integrated profiles 

We first compare the velocity distributions obtained with all the 
data points in the simulation with those obtained by selecting 
only 3% (or 4 x 10 6 ) of these data points. These subsamples 
are of three kinds: a subset of randomly selected positions in 
space and the upper tails of the distributions of the negative 
divergence and vorticity (see Figs. [j]a and |l]b). In the following, 
for the sake of simplicity, we will call them shocks and vortices. 

The choice of 3% of the points is a compromise between 
the value of « 1% dictated by the chemical models and their 
confrontation to the observations and a number large enough 
so that the statistical analysis described below be meaningful. 
The actual values of the thresholds in the vorticity or in the di- 
vergence correspond to this choice. As said above, many of the 
molecular species observed in the cold diffuse medium, require 
large temperatures for their formation. Below gas temperatures 
of the order of 10 3 K, there is an exponential cut-off of the 
speed or the probability of these chemical reactions. Since the 
heating sources are either the viscous dissipation enhanced in 
the layers of large velocity shear at the boundary of coherent 



vortices (Joulain et al., 1998), or the ion-neutral streaming in 



the magnetic precursor of the shocks (Flower and Pineau des 
Forts, 1995), the thresholds in vorticity and in negative diver- 
gence correspond to temperature thresholds. It is only above 
such thresholds in negative divergence (shock velocity) or vor- 
ticity (shear), that the hot chemistry required by the observa- 
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Fig. 2. Velocity distributions obtained with the full velocity 
field and all the points in the simulation (solid), with 3% of 
the data points randomly selected in space (dashed), with the 
3% of the points which have the largest vorticity (dotted) and 
with the 3% of the points which have the largest negative di- 
vergence (dot-dashed). The intensity scale corresponds to the 
number of points in each of the 218 velocity bins. The velocity 
scale is the bin number. 



tions can be triggered in the cold diffuse medium. As seen on 
Fig. [j], the selected thresholds for the vorticity and the diver- 
gence fall in the non-Maxwellian (non-Gaussian) tails of each 
distribution. The sets of structures that we describe therefore 
belong to the non-Maxwellian (non-Gaussian) tails of the cor- 
responding distributions. 

The velocity distributions of each subset (normalized to 
their peak value) is shown in Fig. ^| together with that of the 
whole cube. They are remarkably similar. The only subset 
which provides a slightly (10%) broader spectrum is that built 
on the shocks. A tracer passively advected in a turbulent flow 
(subset of spatially random positions) would therefore carry the 
characteristics of the turbulent field of the bulk of the volume, 
even though its volume filling factor is as small as a few %. 

2.2. Velocity widths of pencil beam and large beam line 
profiles. 

We now compare the velocity width of synthetic spectra ob- 
tained with a large beam (i.e. profiles observed in emission) to 
those obtained in a pencil beam (i.e. profiles observed in ab- 
sorption). To do so we have computed the velocity distribution 
in 4 x 4 x 512 subsamples of the whole cube and of the two 
subsets of largest vorticity and largest divergence. Under the 
approximation that the line radiation is optically thin (Falgar- 
one et al., 1994), we consider that these velocity distributions 
capture the main characteristics of the line profiles. The veloc- 
ity distributions will therefore be called profiles (or spectra) for 
simplicity in what follows. A total number of 1.6 x 10 4 in- 
dividual spectra are therefore obtained across the face of the 
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Fig. 3. PDFs of the velocity widths of the 1.6 x 10 4 individual 
spectra, each sampling 4x4x512 pixels, and therefore mim- 
icking spectra obtained with a beam 128 times narrower than 
that obtained by integrating the emission of the whole cube. 
The spectra are built on (a) the complete velocity field, and on 
only 3% of the data points selected (b) as the most intense vor- 
tices and (c) as the most intense shocks. The vertical bar, in 
each panel indicates the width, at the same level, of the inte- 
grated spectrum over the whole cube. 



cube, for each set. To achieve the comparison of these indi- 
vidual spectra (surrogates of pencil beam spectra) with the to- 
tal spectrum (surrogate of the large beam spectrum), we use a 
width computed at levels 1%, 8%, 16% and 32% of the local 
peak of the velocity distribution (spectrum). Radiative transfer 
affects more severely the velocity domains where the crowd- 
ing is the largest (the peaks of the velocity distributions) and 
this is the reason why we analyze the velocity coverage of each 
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Table 1. Comparison of integrated and individual spectra. 
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level 


Integrated spectrum 
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0.44 
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distribution at levels low enough for radiation to be reasonably 
assumed optically thin. 

We have computed the probability distribution functions 
(PDF) of the widths of the individual profiles at the four levels 
quoted above. They are shown on Fig. || at the 8% level for the 
whole velocity field and the large vorticity and divergence sub- 
sets. The dotted curve shows the Gaussian distributions with 
the same mean and dispersion. On each panel, the vertical bar 
indicates the width of the large beam profile at the same level 
(Autot m Table [l]). The width (Awpop) and peak (i.e. most 
probable value, Av mp ) of the PDFs are given in Table [l] for 
the 8% and 16% levels and differ by less than 20% with the 
level selected. The volume filling factors / v and represent 
the fraction of points included in each sample. 

In all cases, the pencil beam spectra are narrower than the 
large beam spectrum but by only 25 to 45% (column 7 of Ta- 
ble [l]) while the projected sizes of the sampled volumes dif- 
fer by a factor 128. This is due to the fact that although the 
projected size of the volumes sampled by the pencil beams are 
small, the depth along the line of sight has remained unchanged 
(512). Along one dimension at least, the pencil beam samples 
the velocity over large scales. We have made similar computa- 
tions for 8 x 8 x 512 and 32 x 32 x 512 pencil beams and the 
widths of the corresponding velocity distributions remain al- 
most the same. For this reason we believe that the results for a 
pencil beam as small as those of actual observations compared 
to the large beam (i.e. 60"/0.001"=6 x 10 4 , a ratio which cannot 
be achieved by any direct numerical simulation) would not be 
significatively different. As long as the depth of the medium 
sampled by the line of sight is large, the velocity signature 
of the pencil beam profile remains close to that of these large 
scales. 

Then it is interesting to note that the differences between 
the subsets of strong shocks and intense vortices are not 
marked. On average, vortices produce pencil beam profiles 
closer to the large beam profile than do the shocks, but the ef- 
fect in the present simulation is small. 



2.3. PDFs of line centroid increments 

We have computed the line centroid of each spectrum across 
the face of the cube according to the method described in Lis 



et al. (1996). Fig. 2.3 displays the probability distribution func- 



tions of the transverse increments of these centroids for several 
lags. These PDFs are normalized to the rms dispersion of the 
increments. The c omp arison of the PDFs obtained for the full 
velocity field (Fig. |2.3|a) w ith those obtained with the 3% most 



intense vortices (Fig. |2.3|b) and shocks (Fig. |2.3p ) shows that, 
at the smallest lag, the non-Gaussian tails of the PDFs extend 
much further for the latter subsets that for the former. Non- 
Gaussian tails disappear at A = 9 for the full velocity field 
while they persist up to A = 15 for the regions of large diver- 
gence and large vorticity. The non-Gaussian tails of the PDFs 
of line centroid increments have been shown to be associated 
with regions of enhanced vorticity in turbulence (Lis et al., 
1996). This result shows that both the subsets of shocks and 
vortices exhibit pronounced non-Gaussian behaviour in those 
PDFs. 



3. Spatial distribution of the high vorticity/divergence 
subsets 

3.1. Channel maps 

The three channel maps of Fig. || display the surface filling 
factor of the data points in five adjacent velocity intervals for 
the whole cube and the two subsets of vortices and shocks. 
To optimize the visibility of the weakly populated regions, the 
color scale does not span the whole dynamic range and is not 
the same in each velocity interval. At almost all velocities, the 
surface filling factors reach larger values for the shocks than 
for the vortices (see caption of Fig. |]). Thus the shocks are 
more clustered in space than the vortices. Shocks also appear 
as thicker (or more extended) projected structures than vortices. 
Both subsets have a much more pronounced filamentary tex- 
ture than the whole velocity field. Many small scale filaments 
are visible with lengths reaching a significant fraction of the in- 
tegral scale. This result has been already recogniz ed for the re - 



gions of large vorticity by Vince nt & Meneguzzi ( |1991| ; |1 994| ) 



or She, Jackson & Orszag (|1990j) in incompressible turbulence 
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Fig. 4. PDFs of the line centroid in- 
crements for different lags A and (a) 
the full velocity field, (b) the subset of 
most intense vortices and (c) the sub- 
set of positions of largest negative di- 
vergence. 
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Fig. 5. Channel maps of the surface filling factor on each line of sight in (a) the complete velocity field, (b)the subset of the most 
intense vortices and (c) the subset of the largest negative divergence. The color scales are different for each velocity range and 
do not cover the whole dynamic range for the vortices and the shocks. From top to bottom panel, the values of the actual maxima 
are (52, 326, 512, 457, 135) for the complete velocity field, (23, 46, 90, 79, 24) for the vortices and (36, 83, 1 16, 73, 30) for the 
shocks. 
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Fig. 6. Maps of the surface filling factors of (a) the most intense vortices (b) the regions of largest negative divergence. Maps of 
line-of-sight average of (c) the vorticity modulus and (d) the negative divergence. 



and by Porter, Woodward & Pouquet (1998) in compressible 
turbulence. 

For the bulk of the structures, there is no coincidence be- 
tween the patterns delineated by the shocks and the vortices. 
There are exceptions, though, like the structures seen in the 
extreme velocity channels. The small scale patterns delineated 
in these channels by the regions of large vorticity or large di- 
vergence closely follow each other, although the details do not 
exactly coincide. It is remarkable that all the small scale pat- 
terns seen in the extreme velocity channels in the maps of large 
vorticity and divergence are also those of the full velocity field 
(Fig. |]a): in the first channel (top left), the data points of the 
full velocity field are either shocks or vortices or both. 

These maps also show that shocks and vortices are not ran- 
domly distributed but are clustered in space and in velocity. 
There are large voids with only a few vortex or shocks on the 



line of sight. The contrasts are large as indicated by the peak 
values given in the caption of Fig. |[ 



3.2. Integrated maps and surface filling factors 

We now turn to the spatial distribution across the face of the 
cube of the two subsets of shocks and vortices. The set of ran- 
domly selected positions does not have any significant spatial 
structure. It is therefore not shown. The spatial distribution of 
the surface filling factor / s is shown in Fig. ^ for the vortices 
and the shocks. 

Fig. |^ shows that the contrasts are large (w 100 : 1) in 
the two subsets for which the average surface filling factor per 
pixel is only 15.3 (i.e. 3% of the data points over 512 2 pixels). 
There are only 5% and 2% empty lines of sight for the vortices 
and the shocks respectively, i.e. the surface filling factors of the 



8 



Pety & Falgarone: Structure of diffuse molecular gas 



10 



10 



10 



* ID"* 



10 



10 



— high vorticity 

high negative divergence 

random choice 



50 



100 



Fig. 7. Fractional area covered by the structures with a surface 
filling factor above a threshold / sU in the high vorticity (solid) 
and high negative divergence (dashed) subsets. The same curve 
for the randomly selected positions is the dotted curve. 



regions of large vorticity or large divergence are almost every- 
where larger than unity, although they fill only 3% of the whole 
volume. These numbers depend on the thresholds selected for 
the vorticity and divergence. 

The maps of the filling factors (Figs, ^a and |]b) are very 
similar to the maps of the vorticity and of the negative diver- 
gence (Figs. ||c and ^|d) computed with the whole data cube. In 
particular the maxima in vorticity and in negative divergence 
are those of the surface filling factors of the upper tails of the 
corresponding distribution. It means than for the vortices and 
for the shocks, the origin of the maxima observed in projection 
is the crowding of rare events populating the upper tails of the 
distributions. 

We have computed the fractional area covered by structures 
with surface filling factors larger than a threshold / sU . The de- 
pendence of this fractional area on / su is shown in Fig. ^|. Un- 
like randomly selected positions (dotted curve), the high vor- 
ticity/divergence distributions exhibit structures with large fill- 
ing factors which cover only a tiny fraction of the total area. 
The large crowdings of shocks are more numerous than those 
of vortices, but the bulk of the ensemble of vortices (those 
which fill most of the surface, f S Q < 50) cover a fractional 
area slightly larger that the bulk of the shocks. 

4. Comparison with observations 

Statistics on molecular absorption lines in the cold diffuse 
medium are still sparse. We have used the data set of Liszt 



& Lucas (1998) to compare the velocity widths of 9 pairs of 
13 CO(1-0) and 13 CO(2-l) absorption and emission lines to- 
ward extragalactic sources. In 5 cases (B1730, B2013, B2200 
and B2251), the 12 CO lines are simple enough that a compar- 
ison between the emission and absorption profiles is meaning- 
ful. Out of the 14 pairs, two exhibit emission and absorption 



profiles of equal widths, three have absorption profiles larger 
than the emission profiles (this is a possible configuration ac- 
cording to the simulations, see Fig. |]) and nine have absorp- 
tion narrower than emission profiles. The differences are not 
large and the average value of the 9 relative width differences 
(Aw cm — Avabs/AiJem) is only 0.34. This number is very close 
to the values listed in Column 7 of Table [j] for the subset of 
large vortices (although the difference with that of large diver- 
gence is not significant). 

Another interesting property of the molecular line absorp- 
tion measurements is the large scatter of column densities de- 
tected for a given amount of gas sampled by the line of sight, 
iVtot- On average, the observed column density of CH + for in- 
stance increas es alm ost linearly with N to t (see references in 
Joulain et al. ( 1998 )), but the scatter of the values at a given 
value of Ntot is large. Scatters of abo ut a fa ctor 10 are f ound i n 
merging the samples of Crane et al. ( 1995 ) and Gredel ( 1997 ). 
It means that the spatial distribution of the structures bearing 
these molecules is highly heterogeneous, although on average, 
the larger the total column density of gas sampled, the larger 
the column density of molecules. This characteristic may be 
brought together with the fact that the surface filling factor of 
the most intense vortices and shocks have large fluctuations (a 
factor ~ 10) about their average value (15) (Figs. ||a and||b). 

5. Conclusion 



The regions of largest vorticity or largest divergence in com- 
pressible turbulence are small subsets of a whole turbulent ve- 
locity field but we have shown that, despite their small volume 
filling factor (here f v = 0.03): 

(i) they sample the whole velocity field, 

(ii) pencil beam samplings across these subsets have ve- 
locity coverages almost as broad as those obtained with large 
beams. The values (35% to 45%) and the signs of these differ- 
ences are consistent with the observations. 

(iii) the PDFs of centroid velocity increments built on these 
subsets have more extended non-Gaussian wings than those of 
the full velocity field. Their intermittent characteristics subsist 
at larger lags than for the full velocity field. 

(iv) copious small scale structure with large contrasts is 
seen in the maps of the surface filling factor of the regions of 
large vorticity and large divergence. These contrasts are of the 
same order of magnitude as those observed in the column den- 
sities of CH + for instance for a given column density of sam- 
pled gas. 

(v) the surface filling factor of the subsets of high vortic- 
ity/divergence are almost everywhere larger than unity even 
though their volume filling factor is as small as 3% in the 
subsets studied here. Vortices are slightly more efficient than 
shocks at covering the sky: they tend to be more numerous, and 
to form smaller structures which are less clustered in space. 

In summary our study confirms that mild shocks as well as 
intense vortices could be the subsets of the cold diffuse medium 
enriched in molecules and responsible for the molecular ab- 
sorption lines detected in the direction of extragalactic sources. 
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Numerical simulations of hydrodynamical turbulence are 
not ideally suited to test the properties of such subsets, but 
one-fluid simulations of MHD turbulence are not ideal either 
because of the importance of the ion-neutral streaming in the 
triggering of hot chemistry, whether in MHD shocks or in mag- 
netized vortices. More detailed predictions of the phase-space 
structures of the subsets of turbulence where hot chemistry is 
activated in the cold diffuse medium require calculations of 
MHD compressible turbulence which take ion-neutral drift into 
account. 
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